Snakes and ladders: localized solutions of plane Couette flow 
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We demonstrate the existence of a large number of exact solutions of plane Couette flow, which 
share the topology of known periodic solutions but are localized in space. Solutions of different 
size are organized in a snakes- and-ladders structure strikingly similar to that observed for simpler 
pattern-forming PDE systems. These new solutions are a step towards extending the dynamical 
systems view of transitional turbulence to spatially extended flows. 



The discovery of exact equilibrium and traveling-wave 
solutions to the full nonlinear Navier-Stokes equations 
has resulted in much recent progress in understanding 
the dynamics of linearly stable shear flows such as pipe, 
channel and plane Couette flow [H EJ |3j H] . These exact 
solutions, together with their entangled stable and unsta- 
ble manifolds, form a dynamical network that supports 
chaotic dynamics, so that turbulence can be understood 
as a walk among unstable solutions Moreover, spe- 

cific exact solutions are found to be edge states [7], that 
is, solutions with codimension-1 stable manifolds that lo- 
cally form the stability boundary between laminar and 
turbulent dynamics. Thus, exact solutions play a key 
role both in supporting turbulence and in guiding tran- 
sition. 

This emerging dynamical systems viewpoint does not 
yet capture the full spatio-temporal dynamics of tur- 
bulent flows. One major limitation is that exact solu- 
tions have mostly been studied in small computational 
domains with periodic boundary conditions. The small 
periodic solutions cannot capture the localized structures 
typically observed in spatially extended flows. For exam- 
ple, pipe flows exhibit localized turbulent puffs. Sim- 
ilarly, in plane Couette flow (PCF), the flow between 
two parallel walls moving in opposite directions, local- 
ized perturbations trigger turbulent spots which then in- 
vade the surrounding laminar flow [8l [9]. Even more 
regular long-wavelength spatial patterns such as turbu- 
lent stripes have been observed |10j . The known periodic 
exact solutions cannot capture this rich spatial struc- 
ture, but they do suggest that localized solutions might be 
key in understanding the dynamics of spatially extended 
flows. 

Spatially localized states are common in a variety of 
driven dissipative systems. These are often found in a 
parameter regime of bistability (or at least coexistence) 
between a spatially uniform state and a spatially periodic 
pattern, such as occurs in a subcritical pattern forming 
instability. The localized state then resembles a slug of 
the pattern embedded in the uniform background. An 
early explanation of such states is due to Pomeau [IT] , 
who argued that a front between a spatially uniform 



and spatially periodic state, which might otherwise be 
expected to drift in time, can be stabilized over a fi- 
nite parameter range by pinning to the spatial phase of 
the pattern. More recently, the details of this localiza- 
tion mechanism have been established for the subcriti- 
cal Swift-Hohcnbcrg equation (SHE) through a theory of 
spatial dynamics |12l 1131 114] . In one spatial dimension 
the time-independent version of this PDE can be treated 
as a dynamical system in space, in which stationary pro- 
files are seen as trajectories in the spatial coordinate. 
Then localized states correspond to homoclinic orbits to 
a fixed point that visit the neighborhood of a periodic 
orbit representing the pattern. The SHE is equivariant 
under spatial reflections, so the corresponding spatial dy- 
namical system is reversible. There exists an infinite mul- 
tiplicity of reversible homoclinic orbits (i.e., symmetric 
localized states) organized in a pair of solution branches 
which undergo homoclinic snaking. In a bifurcation di- 
agram the two branches intertwine, oscillating back and 
forth within a parameter regime called the snaking or 
pinning region. These are connected by branches of non- 
symmetric states called rungs. Together they form the 
snakes- and-ladders structure of localized states. 

The theory of spatial dynamics also applies to other 
equations in one spatial dimension |15j . but there is no 
obvious extension to higher dimensional PDEs. Nev- 
ertheless, there are remarkable similarities between lo- 
calized states in the simple one-dimensional SHE and 
in other more realistic (and complicated) PDEs. In 
fluid dynamics, homoclinic snaking occurs in driven two- 
dimensional systems such as binary fluid convection |16j 
and natural doubly diffusive convection |17j . Localized 
solutions in these systems exhibit snaking in bifurcation 
diagrams and are homoclinic in that they transition along 
one of the spatial coordinates from a uniform state, to a 
periodic pattern, and back to the uniform state. In shear 
flows, homoclinic snaking has never been observed but 
its existence has been speculated [18]. This speculation 
is supported by the recent discovery of two localized exact 
solutions in PCF by Schneider, Marine and Eckhardt [T5] 
which qualitatively resemble localized states in the SHE. 

The aim of this letter is to elucidate the origin of these 



2 















































































•!■ i 














! 










.... T.tl, 










.. 'I 1 , It. 











(b) 



(c) 



(d) 



2.2 





Mil.' 


■ ,i 




• It.' 






.It, 








. 1 ■ 






; I- 1 




,1 


I mm 






1 mm. 






'.ill 




iil! 


| ] 






',it, " 















FIG. 1: (color online). Localized traveling wave utw (a,b) 
and equilibrium ueq (c,d) solutions of plane Couette flow 
at Re = 400, from 119 . The velocity fields are shown in 
the y — midplane in (a,c), with arrows indicating in- plane 
velocity and the color scale indicating streamwise velocity u: 
blue/green/red correspond to u — —1, 0, +1. The x-averaged 
streamwise velocity is shown in (b,d), with y expanded by a 
factor of three. 
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FIG. 2: (color online). Snaking of the localized utw, ueq 
solutions of plane Couette flow in (Re, D) plane. The spa- 
tially periodic Nagata solution up is shown as well; the utw 
solution connects with it near (131, 1.75). Velocity fields of 
the localized solutions at the saddle-node bifurcations labeled 
a,b,c,d are shown in Fig. [3] The rung branches are shown 
with solid lines connecting the ueq and utw in the snaking 
region; velocity fields for the points marked a, j3, 7 are shown 
in Fig. [4] Open dots on the utw traveling wave branch mark 
points at which the wavespeed passes through zero. 



localized solutions in PCF. We show that the Navier- 
Stokes equations in this geometry indeed exhibit homo- 
clinic snaking, giving rise to localized counterparts of 
well-known spatially periodic equilibria. 

In PCF the velocity field u(x, t) = [u,v,w](x,y, z,t) 
evolves under the incompressible Navier- Stokes equa- 
tions, 

dw. 1 

^ + u-Vu = -Vp+-V 2 u, V-u = 0, (1) 
at Re 

in the domain Q = L X X Ly X L z where x, y, z are 
the streamwise, wall-normal, and spanwise directions, re- 
spectively. The boundary conditions are periodic in x 
and z and no-slip at the walls, u(y = ±1) = ±x. The 
Reynolds number is Re = Uh/v, where U is half the rela- 
tive velocity of the walls, h half the wall separation, and 
v the kinematic viscosity. We treat Re as the control 
parameter and use as a solution measure the dissipation 
rate D — (L x L y L z )~ l / n (|V x u| 2 ) dfl. The laminar pro- 
file has D = 1 while solutions such as those shown in 
Fig. [TJ have D > 1. 

Figure [I] shows two exact solutions of (JT|) at Re = 
400 and f2 = 4tt x 2 x 16-7T, originally identified in [T5] 
for = 4ir x 2 x 8ir. The solutions are localized in 
the spanwise z direction and consist of two to three 
prominent pairs of alternating wavy roll-streak struc- 
tures embedded in a laminar background flow. Fig- 
ure [lja,b) is a traveling-wave solution utw of 
satisfying [tt, v, w)(x, y, z, t) — [it, v, w](x — c x t, y, z, 0), 



where c x = 0.028 is the streamwise wavespeed. Fig- 
ure [l|c,d) is a stationary, time-independent solution 
ueq. The equilibrium ueq is symmetric under in- 
version [u, v, w)(x, y, z, t) = [—u,—v,—w](—x,—y,—z,t), 
and the traveling wave utw ua s a shift-reflect symmetry, 
[u, v, w] (x, y, z, t) — [u, v, — w]{x + L x /2, y, — z, t). These 
symmetries ensure that neither ueq nor Utw drifts in 
the localization direction z. 

To continue these solutions in Re, we combine a 
Newton-Krylov hookstep algorithm [2U] with quadratic 
extrapolation in pseudo-arclength along the solution 
branch. The Navier-Stokes equations are discretized 
with a Fourier-Chebyshev-tau scheme in primitive vari- 
ables and 3rd-order semi-implicit backwards differenti- 
ation time stepping. Bifurcations along the solution 
branches are characterized by linearized eigenvalues com- 
puted with Arnoldi iteration. The computations were 
performed with 32x33x256 collocation points and 2/3- 
style dealiasing and checked against computations with 
(3/2) 3 more gridpoints. The numerical software is avail- 
able atlwww7channelflow.org 

The bifurcation diagram in Fig. [2] shows the Utw 
and Ueq solutions from Fig. [TJ under continuation in 
Reynolds number. As Re decreases below 180, the so- 
lution branches snake upwards in dissipation D; that is, 
they pass through a sequence of sub- and supercritical 
saddle-node bifurcations which nearly line up, creating 
a large multiplicity of localized solutions in 169 < Re < 
177. Each saddle-node bifurcation adds structure at the 
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FIG. 3: (color online). Localized traveling wave utw (left) and equilibrium ueq (right) solutions of plane Couette flow at 
points marked on the solution branches in Fig. [5] (a,c) show the velocity fields of utw at its first and second saddle-node 
bifurcations, moving up each branch from lower to higher dissipation D; similarly (b,d) for ueq. (e,f) show the a;-averaged 
velocity, with in-plane velocity indicated by arrows and streamwise velocity by the colormap. The marginal eigenfunctions at 
the saddle- node bifurcations (c,d) are shown in (g,h). 
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FIG. 4: (color online). Localized solutions of plane Cou- 
ette flow along a rung branch, for the points marked a, /3, 7 
in Fig. |2j and plotted in terms of x-averaged streamwise 
velocity (u) x (y,z) as in Fig. ljb,d). (a) shows the begin- 
ning of the rung solution on the ueq branch with symmetry 
(u) x (y,z) = -{u) x (-y,—z). Midway along the rung, (/3) 
is non-symmetric. The rung terminates at (7) on the utw 
branch with symmetry (u) x (y,z) = {u) x (y,— z). 



edges (fronts) of the localized solution while preserving 
its symmetry. This spatial growth is illustrated in Fig. 
[3] which shows the velocity fields at several points along 
the snaking branches. For example, Fig. [3]ja) shows the 
y = midplane of utw &t the saddle-node bifurcation 
marked a in Fig. [2] Continuing up the solution branch 
from (a) to (c), the solution gains a pair of streaks at 
the fronts while the interior structure stays nearly con- 
stant. The marginal eigenfunction associated with the 



saddle-node bifurcation at (c) is shown in Fig. |3jg); it 
is weighted most heavily at the fronts of the localized 
solution and has the same symmetry, so that adding a 
small component of the eigenfunction strengthens and 
slightly widens the fronts, whereas subtraction weakens 
and shrinks them. 

The spanwise wavelength of the interior structure of 
the localized solutions is approximately £ z = 6.8 ~ 2n. 
This value is selected by the fronts that connect the in- 
terior streaks to the laminar background, and it does not 
seem to vary across the snaking region or when compared 
between the two branches. The streamwise wavespeed 
c x of the traveling wave solution varies along the branch 
and in fact changes sign several times. Points at which 
c x = are marked in Fig. [2] with open circles. The point 
marked a has c x — 0.0062. Rotation about the z axis 
generates symmetric partners for both Uxw an d ueq; 
for the former this results in a streamwise drift in the op- 
posite direction. Thus the utw and ueq curves in Fig. 
[2] represent four solution branches. The lower branches 
of both can be continued upwards in Reynolds number 
past Re = 1000. 

Figure [2] also shows six rungs of non-symmetric ex- 
act localized solutions. These bifurcate from the snaking 
branches close to the saddle nodes and are associated 
with marginal eigenfunctions whose symmetry does not 
match the base state. Each rung connects to both the 
Ueq and the u TW branch so solutions along the rungs 
smoothly interpolate between the two symmetry sub- 
spaces, as illustrated in Fig. [4] The rung solutions travel 
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in z as well as x, but with z wavespeed three orders of 
magnitude smaller than c x . 

Due to the finite extent of the domain, the struc- 
tures cannot grow indefinitely, and the snaking behav- 
ior must terminate. As in other problems of this type, 
the details of this termination depend on a commensu- 
rability condition between the spanwise wavelength of 
the streaks within the localized solutions and the span- 
wise domain |21j . At L z = 167r, Uxw connects at 
(Re,D) = (131,1.75) to the spatially periodic Nagata 
equilibrium with wavelength t z = 2ir . The ueq branch 
does not appear to connect to any periodic solution. In- 
stead, when this branch exits the snaking region its ve- 
locity field contains a localized defect that persists under 
continuation up to at least Re = 300. At other values 
of L z , both branches might cither not connect to a peri- 
odic solution or connect to a different periodic solution. 
For example, in Ref. [19] it was shown that at L z = 8ir 
both the utw an d ueq branches terminate on a branch 
of spatially periodic solutions, though that choice of L z 
was too narrow to allow the snaking structure to develop. 

We have shown that homoclinic snaking in wide plane 
Couette channels gives rise to a family of exact local- 
ized solutions with internal structure similar to the pe- 
riodic Nagata equilibrium. Thus, as recently specu- 
lated [TBI EH] > the localization mechanism studied in the 
SHE carries over to linearly stable shear flows, where 
it cannot be associated with an instability of the uni- 
form state [22 . Physically, the localized states studied 
above consist of fronts pinned to the periodic Nagata 
equilibrium. The periodic structure is formed by pairs of 
counter-rotating, streamwise-oriented roll-streak struc- 
tures, which also characterize other exact solutions linked 
to transitional turbulence in small domains. Therefore, 
localized versions of the other known exact solutions 
should also exist and, together with their heteroclinic 
connections, support localized turbulence. In this sense 
the localized solutions studied here are a first step to- 
wards generalizing the dynamical systems picture for tur- 
bulence to extended flows. 

Turbulent spots and stripes that are tilted against the 
flow direction suggest the existence of exact solutions lo- 
calized in both spanwise and streamwise direction. Al- 
though a theory for localization in two spatial dimen- 
sions is not yet available, numerical studies of the SHE 
show that snaking does carry over to solutions localized 
in two dimensions [23) . In PCF it is however not known 
if the same mechanism also generates fully localized ex- 
act solutions, because a bifurcation analysis would first 
require a fully localized solution to start the continua- 
tion. Such a solution is unfortunately not yet available. 
Nevertheless edge calculations both in pipe flow [211 12"5] 
and in extended plane Couette cells [JjJJ [5SJ [57] yield 
localized structures that show very mild dynamic fluctu- 



ations, which suggests the existence of simple underlying 
fully localized exact solutions. 
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